Landau-Zener Tunneling of Bose-Einstein Condensates in an Optical Lattice 
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A theory of the non- symmetric Landau-Zener tunneling of Bose-Einstein condensates in deep 
optical lattices is presented. It is shown that periodic exchange of matter between the bands 
is described by a set of linearly coupled nonlinear Schrodinger equations. The key role of the 
modulational instability in rendering the inter-band transitions irreversible is highlighted. 



Introduction. After the first realization of Bose- 
Einstein condensation (BEC) in an optical lattice 
(OL) ^J, and subsequent experimental studies of Bloch 
oscillations and Landau-Zener (LZ) tunneling of mat- 
ter [2[, numerous theoretical efforts were made to un- 
derstand and describe these phenomena. In particular, 
the effect of the nonlinearity on LZ tunneling has been 
recently studied in Refs. 

SB- 

Those studies, however, 
are essentially based on the assumption of smallness of 
the lattice potential in comparison with the recoil en- 
ergy, which corresponds to a narrow gap in the spectrum 
of the underlying linear system - i.e., to the case where 
one can employ a two-mode approximation, where the 
modes are plane waves. Meanwhile, recent direct exper- 
imental observation of the LZ tunneling of a BEC in an 
OL 0, E| was performed for the range of the potential 
amplitudes extended till several recoil energies (a deep 
OL), challenging the existing theory 7]. In the latter 
case, the fine structure of Bloch states becomes of crucial 
importance and leads to a model different from the one 
proposed in 0, 3 . Although the basic principles of the 
LZ tunneling are hosted in the linear band structure, the 
nonlinear interaction makes the phenomenon very differ- 
ent from the one known in solid state physics. Indeed, in 
the presence of nonlinearity, the Bloch states can become 
unstable |8| and localized states can appear in the gaps 
between bands (gap-solitons) which can dramatically af- 
fect the atomic transfer rate. 

The aim of the present paper is to provide a theoretical 
understanding of the nonlinear interaction and LZ tun- 
neling between the lowest bands of a periodic potential. 
This is particularly relevant for a deep OL, where the 
two edges of the lowest forbidden zone have essentially 
different properties. This approach reveals not only the 
role of modulational instability (MI) on the atom trans- 
fer between zones, but also the effect of (partial) sup- 
pression of the MI by the lattice acceleration. Such an 
understanding is important for appreciating the ramifi- 
cations of this key linear phenomenon in the presence of 
nonlinearity and its connection with critical nonlinear in- 
stabilities such as the MI. It is also particularly relevant 
for interpreting recent experiments in this system [3, |£| . 



Theoretical Approach. We start with the Gross- 
Pitaevskii equation for the macroscopic wave function 

= -—A^ + V t rap(r)^ + Vi aU (r)^ + go\^\ 2 ^, (1) 

where r = (x,r±), go = 47rh 2 a s /m, a s is the s-wave 
scattering length, and m is the atomic mass. The trap 
potentials are given by V tr ap{ v ) — ^ (^^ 2 + ^o 7 "!) an d 
by Viatt(r) = cos 2 [x - § t 2 )] where u x and co are 
the harmonic oscillator frequencies, A is the laser wave- 
length and a is the acceleration of the OL. 

In the following we restrict to the weakly nonlinear 
regime which can be quantified as e = (2y/2^k)~ 1 <C 
1 where £ = (87rn|a s |) -1 / 2 is the healing length, k = 
2^, and n is the mean atomic density. To this end, we 



introduce new variables X 
and T 
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is the recoil energy, 
the dimensionless parameters v = le j^ H •> 7 = 7%w£& -> anc ^ 



h z k z 

2m 



v a = a =sign(a s ), and a wave function \I/(T, R): 
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Then, Eq. Q 

can be rewritten in the form 



id T V = - (C x + £±)* + jX (l + -fT 2 ) * 

v 2 



(3) 
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where Cx = —d x + vcos(X) and Cj_ = — 

For weak nonlinearity the energy of the two-body in- 
teraction is much less than the recoil energy and Eq. (|3J) 
can be reduced to a system of coupled one-dimensional 
nonlinear Schrodinger equations. To show this we intro- 
duce a set of variables t\ = e z T, x\ = e z X, (Z = 0, 1, ...) 
which are considered independent, and look for a so- 
lution in the form \I/ = ipi + eip2 + where ipj = 
ipj(to, xo, R_l). After substituting these expansions 
into (PJ, we single out equations of the same order in e as: 



[id t0 -%(C X0 +C±)]il; j = Fj, where F x = 0, F 2 
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etc., and C XQ = —d 2 Q + vcos(xq). We then consider the 
eigenvalue problems C Xo (f aq (x ) = S a (q)(p aq (x ), and 
>Cj_C/x_l = l J '±(n±i where a characterizes the zone, q is 
a wave- vector in the first Brillouin zone (BZ), and fi± 
refers to the transverse quantum numbers. 

Focusing on the tunneling (i.e., particle transfer) be- 
tween the two lowest bands, straightforward algebra 
shows that if all atoms are initially concentrated in the 
lowest transverse state, the lattice acceleration does not 
cause (to leading order) transitions among the transverse 
levels. This allows us to use (± = (^) 1/2 exp(-fR 2 ± ) 
and look for the wave function in the form 

fa = e'^^aiAr { Xl M)vi + A 2 (x u t 2 )cp 2 ] (4) 

where we denoted with cpi y2 the modes bordering the BZ 
(i.e. the ones with q = tt/2). Further simplification can 
be achieved by assuming that the band structure is well 
pronounced, i.e., that the condensate is wide enough in 
the ^-direction, huu x <C Er. This is expressed through 
the rescalings v x = e 2 v and 7 = e 2 7 where v ~ 7 ^ 1, 
allowing one to neglect the term proportional to jv 2 T 2 . 

We now use the multiple scale technique (see e.g., 0) 
where we take into account the existence of two modes 
due to tunneling. This leads to the system 

id t2 A x + (2m 1 y 1 d 2 xl A 1 - PnAx - f3 12 A 2 
HD 2 x 2 /2)A 1 - (xnl^il 2 + Xi2|A 2 | 2 )A! = 0, 
id t2 A 2 + (2m 2 )- 1 d 2 Xi A 2 - p 21 A ± - (3 22 A 2 w 
+ (v 2 x 2 /2)A 2 - (x2i|A!| 2 + X 22\A 2 \ 2 )A 2 = 0. 

Here m" 1 = | d are the inverse effective masses, 

\(pi(x)\ 2 \(pj(x)\ 2 dx (6) 

7 r 

Pij = ^ J x(pi(x)(pj(x)dx. (7) 

Recalling that cpj border the lowest gap, one concludes 
that i) mi < and m 2 > 0; ii) flu = (5 22 = and 
P12 — P21 — A and hi) X12 = X21 = X (below we use the 
notation Xj = Xjj)- Then, we end up with the system 

id t2 A a = -(2m a )- 1 dl 1 A a ^(3A 3 _ a + (9 2 x 2 /2)A a 

+ (Xa\A a \ 2 + X \A 3 -a\ 2 )A cn a =1,2. (8) 

Approximating the Bloch functions by <p a (x) — 

>/2cos( a!+(1 - a),r ), we find that x = Xj = 3%, /3 = 7- 
At this point, it is worth to emphasize the differences 
between our model and the one used in Refs. 0,0. First, 
in a generic case the nonlinear terms responsible for self- 
phase and cross-phase modulation appear with different 
coefficients, respectively xi and x? originating from the 
different structure of the Bloch functions at two edges of 
the gap. Unlike the phenomenological model considered 
in Ref. these terms appear in (J2J) with the same sign. 



Second, our model includes the group velocity disper- 
sions (i.e. the inverse effective masses) in explicit form, 
a feature that does not occur in the approximation of a 
shallow OL. Third, the time does not enter explicitly the 
equations for the envelopes because it is scaled out by 
the ansatz (0 using smallness of 7^, making the result- 
ing system conservative. Finally, our model takes into 
account the trap potential (an essential ingredient of the 
experimental setup jl,0|)- 

The system JJJ leads to a number of conclusions: 

(a) In the presence of nonlinearity the LZ tunneling is 
an asymmetric effect. This stems not only from the fact 
that the coefficients of the two equations are different i.e., 
|mi| _1 7^ m^ 1 and xi 7^ X2, but also form the fact that 
the condensate at one gap edge is modulationally stable 
[edge 2 (1) in the case of a positive (negative) scattering 
length] and modulationally unstable at the other edge || . 

(b) If, however, one considers the evolution of un- 
modulated Bloch waves (i.e., the case d Xl Aj = 0) then 
the inverse effective masses do not appear in the descrip- 
tion of the phenomenon and in a rather common situa- 
tion when xi ~ X2, the LZ tunneling becomes symmetric. 
This case can be explicitly treated analytically. The anal- 
ysis of the resulting equations for the populations of the 
two bands demonstrates the existence of oscillations of 
matter between the bands. However, as this case ignores 
the important feature of MI, we do not discuss it here. 

(c) Although Eqs. © are coupled for 7 = 0, the 
number of particles Mj = / \ Aj\ 2 dxi, is separately pre- 
served (i.e. dMj/dt 2 = if 7 = 0) and LZ tunneling 
cannot occur in this case. On the contrary, for 7 7^ 
one has the conservation of the total number of particles 
J\f = J\f\ + M 2 : dM I dt 2 = and LZ tunneling can occur. 

Numerical Results. We define the rate characterizing 
tunneling as r(t) = Mi/M. In the numerical experi- 
ments, we have simulated the 2-state transition model 
of Eqs. JHJ) with parameters adapted to the experi- 
ments reported in Ref. p. In particular, we consider 
V = 2.2E R w h x 10 3 Hz, corresponding to v = 0.1375; 
hence, m^ 1 ^ — 1.5 and m^ 1 ~ 2.5 (as obtained from the 
band structure of the respective potential) . A weak mag- 
netic trapping for rubidium atoms has been implemented 
by choosing uo x ~ 2tt x 55Hz and uoq w 2tt x 200Hz, which 
corresponds to v x = 0.03 and v$ = 0.1. The nonlinear 
and linear couplings have been obtained from Eqs. (0- 
0. In particular, x = 0.008 and Xj = 0.024, while /?, 
similarly to the setting of 6], was varied (typically) in 
the interval /3 G [0.2,1.4]. Furthermore, the calculations 
were also performed for solutions of various amplitudes 
to emulate the variation of the nonlinearity in p. 

We first examined cases where Ai(x, 0) = and the 
second component is an exact solution of Eq. (0 with 
a = 2 and in the absence of linear coupling. Typical re- 
sults for the transfer rate are shown in Fig. ^ At the 
early stages of the dynamical evolution, one observes pe- 
riodic particle transfer between two bands, reflecting the 
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FIG. 1: The tunneling ratio r(t) — N\/N is given for three 
different amplitudes and /3 = 0.8 on the top three panels: 
max(A2(x, 0)) = 3.22 (first from the top), 4.83 (second) and 
6.46 (third). It is also shown for max(A2(x,0)) = 4.83, for 
three different couplings (3 = 0.2 (fourth panel), 1.4 (fifth) 
and 10 (sixth). 



FIG. 2: The first two panels from the top show the spatio- 
temporal contour plots of |Ai(x,£)| 2 (first from the top) and 
|A2(x, t)| 2 (second panel). Initially (i.e., at t = 0) the second 
component is populated [max(^2(x, 0)) = 4.83], while the 
first is not. The remaining three panels show three snapshots 
at different times of the spatial profile of the first (solid line) 
and of the second (dashed line) component. In the middle 
and bottom panels the development of the MI is evident. 



fact that the probabilities of direct and reciprocal tun- 
neling are the same (an effect not taken into account by 
the simplified qualitative picture in |6| but also predicted 
by the model of Ref. 4]). However, after some critical 
time, t cr , the behavior of the system changes drastically: 
the amplitude of the tunneling rate oscillations decreases 
substantially, indicating that the system never returns 
back to the state where all atoms are condensed in only 
one band. As one can see from the top three panels in 
Fig. EI the critical time is smaller for larger initial atomic 
populations of the second zone i.e., for larger nonlinear- 
ity. The remaining panels of Fig. Q] display two other 
important features: (i) the frequency of oscillations in- 
creases with the lattice acceleration and (ii) the critical 
time also increases with the lattice acceleration. The first 
of these features seems to be intuitively clear: larger lin- 
ear force corresponds to larger probability of tunneling 
and thus to more intensive interchange of atoms between 
the bands. The second item, as well as the existence of t cr 
itself, are less evident. In order to appreciate these fea- 
tures, we computed the evolution of the density profiles 
corresponding to the two bands. Fig. [2]shows the spatio- 
temporal contour plot (as well as snapshots of time) evo- 
lution of the two components for an initial condition with 
max(A2(x,0)) = 4.83. Approximately at t = 35 (for the 
case e 2 ~ 0.1, corresponding in physical units to approx- 
imately 90 fis), MI starts to develop in the system and 



at t = 60 one can already see a well pronounced peri- 
odic structure of the distribution of atoms in the first 
band, which at t = 90 develops into a number of soli- 
tary pulses, as is typical for MI. The emergence of stable 
localized states eventually renders the transfer of mat- 
ter between the two zones irreversible. The time of the 
development of soliton-like structures coincides with the 
critical time observed in Fig. (second top panel). This 
explains the change in the dynamics of the tunneling rate 
reported above. Namely, the particle exchange between 
the two bands is suppressed after solitonic (train) struc- 
tures are developed from the "unstable" band (the lowest 
band in our case). Notice that the stronger the nonlin- 
earity is, the larger is the region of unstable wavenumbers 
and equivalently the shorter is the time for the instabil- 
ity development, thus explaining the differences in the 
tunneling rates shown in Fig. ^ One can also explain 
the dependence of r(t) on the lattice acceleration (linear 
force). Indeed, while in the "unstable" band, the non- 
linearity tends to fragment the condensate, immediately 
after transfer to the "stable" band the nonlinearity leads 
to dispersion of the wave packet. By increasing the lin- 
ear force (/?) one increases the frequency of oscillations 
of r(i), i.e., decreases the time atoms spend in the "un- 
stable" band, thus effectively reducing the MI growth 
rate. This explains the increase of the critical time t cr 
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at which MI develops, observed for increasing (3 in the 
three bottom panels of Fig. ^ notice from the bottom 
panel that for very large j3 the effect saturates. In all the 
cases we observed rather weak dependence of the final 
upper-to-lower tunneling on the nonlinearity which cor- 
roborates the results of Ref. p. Our result of the tun- 
neling rate however is larger than the experimentally ob- 
served, which, apparently, is an artifact of the two-mode 
model; the latter does not account for the atoms tunnel- 
ing to higher bands, thus reducing the effective number 
of atoms participating in inter-band transfer with respect 
to the experiment. 

The case in which all atoms are initially in the low- 
est band can be treated similarly. Since, in this case, 
the homogeneous distributions are unstable, it is natural 
to start with an initial condition which is a gap soliton 
resulting from MI and containing the same number of 
atoms as in the previous case. To understand the phe- 
nomenon, we use as initial condition a soliton which is 
somewhat displaced from the bottom of the potential so 
that the tunneling matter will be spatially separated from 
the pulse (due the dynamics in the trap) and therefore 
easily distinguishable. This setting is illustrated in Fig. 
3. As is shown in the snapshot of the bottom right panel 
of the figure, the wave packets on the left are well sepa- 
rated from the localized pulse. These represent the mat- 
ter that has tunneled to the second band. Our numerical 
results indicate that to each burst of matter in the second 
band, there corresponds one in the first band. Also no- 
tice that the burst of matter in the second band (see the 
top right panel of Fig. 3) is becoming more and more ex- 
tended as one expects for a Bloch wave of the stable band. 
The separation of the matter strongly decreases the inter- 
action between the bands (the corresponding wave pack- 
ets have practically zero overlapping), thus leading to 
the formation of a plateau in the tunneling rate 9]. The 
asymptotic tunneling rate achieved in this case may be 
different from the one obtained in previous case (cf., e.g., 
Fig. Q); this further supports the existence of an essen- 
tial asymmetry induced by the nonlinearity in the LZ 
tunneling |lflj . 

Conclusions. We have provided a systematic deriva- 
tion of a two-band model reduction relevant to the in- 
teraction of the lowest bands of a periodic potential. 
A theme of particular interest, in this context, is the 
Landau- Zener tunneling of a BEC between two bands of 
an optical lattice, a recurring theme of theoretical and ex- 
perimental importance. While a number of our results, 
like the sensitivity of the effect with respect to initial 
distributions, the oscillatory behavior of tunneling rate, 
and the qualitative symmetry with respect to exchange 
of bands 1 <-> 2 with simultaneous change of the sign of 
the scattering length corroborate the predictions of 
we have illustrated important differences from previously 
explored models and have highlighted the asymmetric 
nature of the tunneling and the important role of the 
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FIG. 3: The top left and right panels show space-time contour 
plots of the 1st and 2nd band wavefunctions respectively. The 
bottom left shows the temporal evolution of the tunneling 
ratio r(£), while the bottom right shows a snapshot of the 
solution at t — 50. Solid line: |^4i| 2 ; dashed line: |^b| 2 ; dash- 
dotted: parabolic potential. 

nonlinearity through the modulational instability in the 
process. These results illustrate the importance of (ex- 
perimentally) studying this system for longer times sub- 
stantially larger than the critical time (and/or stronger 
nonlinearities) to appreciate the role of the nonlinearity 
in modifying the classical picture of the Landau-Zener 
tunneling. 
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